Structure, Stresses and Local Dynamics in Glasses 

Tamar Kustanovich and Zeev Olami 

Department of Chemical Physics, 

The Weizmann Institute of Science, Rehovot 76100, Israel 



On 

c 

o 
m 

c 

C/5 



The interrelations between short range structural and elastic aspects in glasses and glass forming 
liquids pose important and yet unresolved questions. In this paper these relations are analyzed 
for mono-atomic glasses and stressed liquids with a short range repulsive-attractive pair potentials. 
Strong variations of the local pressure are found even in a zero temperature glass, whereas the largest 
values of pressure are the same in both glasses and liquids. The coordination number z{ J) and the 
effective first peak radius depend on the local pressures J's. A linear relation was found between 
components of site stress tensor and the local elastic constants. A linear relation was also found 
between the trace of the squares of the local frequencies and the local pressures. Those relations 
hold for glasses at zero temperature and for liquids. We explain this by a relation between the 
structure and the potential terms. A structural similarity between liquids and solids is manifested 
by similar dependencies of the coordination number on the pressures. 

PACS numbers 61.43 
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I. INTRODUCTION 

Short range order in glasses and liquids has a set of dif- 
ferent interrelated aspects: structure, local stresses and 
the elastic properties (e.g. local frequencies). Though a 
lot of attention has been devoted to the structural point 
of view and to frustration effects [1-7], almost no atten- 
tion was given to the local elastic properties l| and to the 
existence of a wide range of local stresses in such systems 
P, E5Q |. In a liquid state, there is a great variety in the 
local pressures felt by atoms, as a result of the variations 
in local environments M. When a glass is created from 
the liquid state, its atoms will reach one of the local min- 
ima of energy (defined by zero force) , but not the global 
energy minima (corresponding to the crystal). There- 
fore, the local pressures in the compressed liquid will not 
be fully relaxed in the quench. Furthermore, the exis- 
tence of different stressed local environments will result 
in varying local elastic constants in compressed liquids 
and glasses. 

Crystals and glasses are notably similar, since both 
have a fixed structure around which the zero point dy- 
namics can be discussed. Yet, there is a basic structural 
difference between them. In crystals, due to the crys- 
talline order, there is a unique lattice with well defined 
distances between sites and periodic long range order 
that corresponds to the structure at minimum energy. 
The dynamical properties and elastic constants are ob- 
tained by expanding the energy around that unique ref- 
erence lattice. In glasses there is only short range order, 
which is determined by the zero force condition and by 
the fluctuations in the local pressures. On larger scales, 
the structure is homogeneous. 

It is important to realize that the knowledge of struc- 



tural aspects of the glass cannot supply any information 
about the glassy dynamics. Furthermore, to understand 
the glassy structure one has to consider information re- 
lated to the internal stresses. Hard sphere glasses can 
supply a useful first order description of structure in 
glasses 0]. However, packings with very similar structure 
can have a very different elastic description, so knowledge 
of the forces is essential. Even if some glassy structure 
is known, a redefinition of the nearest neighbors interac- 
tions can completely modify the spectra. In this paper 
we explore this problem, focusing on the relations be- 
tween the local dynamics, stresses and structure. We 
have found that the fluctuations in structure and in the 
force are combined together to create a linear relation be- 
tween the local elastic constants and the local pressures. 

This kind of analysis can enhance the understanding 
of structure in glasses. It can be used to predict the 
variability of stress and elastic constants in glasses, which 
can then be used to give predictions on the global spectra 
and other glass characteristics. 

In this paper we discuss short range mono-atomic 
glasses (SRMA glasses) and stressed liquids. In section 
III it is shown that local pressures are distributed in a 
wide stable distribution, which is related to the pair dis- 
tribution function. There is a common upper cutoff for 
the local pressure distributions in glasses and liquids. We 
have found that the variations in pressure and structure 
are correlated and that features of the non-stable liquid 
state are sustained in the glass. 

Section IV describes correlations between stresses and 
local dynamics. We have found a linear relation between 
local pressures and the elastic constants. This relation 
can be explained by combining the structural information 
discussed in section III with knowledge of the potential. 



Thus, this linearity is a strong constraint on the structure 
of 6-12 Lennard-Jones glass. Furthermore, it seems to be 
a general feature of SRMA glasses (see section V) . 



II. ELASTIC DESCRIPTION 

First, we define our basic points of interest. Consider 
a solid at equilibrium state at zero temperature. There 
are no forces on atoms and no external stresses on the 
boundaries. In such case, the average stress and, in par- 
ticular, the average pressure are zero. Note that to ensure 
absence of stresses on boundaries, one needs to consider 
only potentials with attracting and repulsive terms (e.g. 
6-12 Lennard-Jones potential) with density adjusted so 
that total pressure equals zero. In crystals, the aver- 
age stress (pressure) and the deviations from it are zero. 
This is not so in amorphous solids: even if the amor- 
phous system under consideration is left with a zero av- 
erage stress, the local stresses can be very different from 
zero @||li]]. 

Consider a discrete system where atoms i and j in- 
teract through a pair potential U(Rij), where R/j is the 
distance between two atoms. The forces between parti- 
cles are defined by: 



.h = J2u'(Rv) 
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(1) 



where U'(Rij) = dUij/dRij and R\j = (Xij,Yij, Zij). In 
a glass at zero temperature: 



ji= J2 u '^) 



(6) 



as the local pressure. The distribution function of the 
local pressures is denoted as P(J). 

The local dynamics of an atom is defined by the matrix 
A % = J2 ^a/3 which describes the second derivatives of 
energy in particle i while other atomic positions are fixed. 
The A y 's are defined as: A l ° a/3 — d a dpU % i , where a 
and j3 denote components of a vector. For radial pair 
potential 
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where U ii — d 2 Uij/dRlf, and I is identity matrix. Note, 
that there are terms in A which are related to the first 



gradients 



The trace of this tensor is: 



K % = trace{Ai) = Y^l U ij 
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A liquid of the same volume has velocities and the 
zero force condition does not hold. However it is still 
interesting to use the previous elastic descriptions. The 
stress, pressure and local elastic frequencies can still be 
measured, even though they do not have such a direct 
relationship to the spectra. Doing so, we ignore the effect 
of velocities on the pressure. 
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The site stress tensor can be defined as fl( 
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The average stress field around atom i can be defined 
as the average stress in a cell of volume Vt 



Oi = — 

Vi 



(4) 



Usually this cell is supposed to be the Voronoi cell, how- 
ever this is not a necessary assumption [flOp]. At equi- 
librium, it follows that 
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all— atoms 



However there is no such constraint on the various Oi. 
The same requirement holds for the pressures. We define 
the term 



III. STRESS AND STRUCTURE 
A. Systems Studied 

We considered systems with 4096 atoms interacting 
though 6-12 Lennard-Jones potential, using the values 
e/fcfc = 120 and a = 0.34nm, which correspond to Argon. 
Periodic boundary conditions were assumed. The density 
p = 1.763gr/cm 3 was chosen to ensure no external pres- 
sures at equilibrium. The samples were prepared from a 
compressed liquid at 120K, by using a molecular dynam- 
ics simulation ]13fl , either by a steepest decent quench 
to zero temperature glass, or by coupling to a heat bath 
at temperatures of 1 and 8K. Details of glass prepara- 
tion by various methods are summarized elsewhere [ p_ 1| . 
All those glasses were quenched to zero temperature, to 
ensure the existence of an elastic description. 

We calculated the local Si, Ji, Aj and K~i (eq. 3, 6, 7 
and 8) on various samples. Some initial calculations were 
made with other potentials as well |L7| (see comments in 
section V). 



B. Local Stresses and Structure 

A fully relaxed glass with zero global pressure has large 
variation of local pressures J's. Therefore, we first look 
at the width of distribution and discuss the essential fea- 
tures of the local pressures in our system. The distribu- 
tions of J for a glass at T=0 and for liquid at T = 120Jf 
are presented in Figure 1. 




FIG. 1. P(J) vs. J for liquid at 120K (x) and for glass at 
OK (■). J's units are e/a, where e and a are Lennard- Jones 
parameters. 

The distribution in glass, P g i(J), is wide, in the sense 
that its width is of the same scale as the typical J's for 
which the system will fracture (at uniform expansion) 
[ JLjJ . The average of J is not zero even though the sys- 
tem is under a zero pressure (see eq. pi) . There is a sharp 
peak in P g i(J) at J = 0. This peak is due to about 1 — 2% 
of particles that were found to be in an 'ideal' equilibrium 
state (they are located in exact potential minima). Those 
minima are due also to neighbors which are far from the 
first peak. The preparation methods do not affect sig- 
nificantly the width of those distributions. In fact P(J) 
is as stable as the radial pair distribution function g(r). 
Changes occur only due to crystallization. 

The liquid distributions Pu q (J) are wider || and the 
zero frequency peak is much smaller. Yet they are similar 
in form to the glassy distributions. A surprising aspect is 
that the highest local pressure J ma x is identical in the 
liquids and the glasses distributions. It is impossible to 
generate a site with a higher local pressure. We discuss 
this point later. The shapes of P g i(J) and Pu q {J) can be 
re-scaled on each other, using the relative width of the 
distributions (or the average pressure) and the maximal 
value J m ax- 

In glasses and liquids, the usual description of the ra- 
dial structure is given through the radial pair distribution 



function g(r), which is defined by the average number of 
neighbors at distance r. This function is not sensitive 
to variations in the local quantities. An integration of 
the type J av = J g{r)U'(r)dr will give only the aver- 
age local pressure in the glass, but will not provide any 
information about the variability of J or of the other pa- 
rameters. Conditional radial pair distribution gj{r) is 
defined as the average number of neighbors at distance 
r for a specific given J. Since g(r) = J P(J)gj(r)dJ, all 
the information about the variability of J is given indi- 
rectly in g(r). However, the structural effects of the local 
pressures can be observed directly, using the conditional 
radial distributions gj(r). 
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FIG. 2. gj(r) vs. r for 3 averaged values of J - solid curve 
for compressed local environments, dotted curve for low pres- 
sures (small values of \J\) and dashed-dotted line for stretched 
local environment (high values of J), r's in units of 2 ' 6 <r. 

The conditional distribution functions gj(r) are shown 
in Figure 2. As expected, J is dominated by contribu- 
tions from the first peak while there is an almost con- 
stant contribution from further neighbors. This is due 
to the fact that the pressures are much smaller at large 
distances and the distribution is smoother. A reasonable 
approximation to J is: 



J J fp T" Jb g 



(9) 



where we denote f p and b g as the first peak and further 
neighbors contributions, respectively. A numerical esti- 
mate for Jb g is (11.5 ± 1.5)e/<7. Though the main effect 
on the local pressure arises from the first peak in gj(r) 
there are large structural effects in the second peak. This 
finding is related to midrange order in glasses, which was 
discussed extensively in the last decade Hjlq] and will 
be addressed in a forthcoming paper [H} . 



The first peak is described by two major structural 
parameters: the coordination number and the average 
peak position, given by the following definitions: 
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where the cutoff of those integrals is defined by the end of 
the first peak. Two major effects arise from the analysis 
of the dependence of z(J) and R(J) av on J: a shift of the 
average position of the peak R(J) av and a change in the 
coordination number z(J). z(J) is presented in Figure 
3. We have found a square dependence of z(J) on J for 
medium and high local pressures. 




FIG. 3. (a) z(J) vs. J for two different samples of glass 
at OK. The plot has a good quadratic fit for J > 25 (except 
for the upper point which is statistically meaningless as it 
represents ~ 0.006% of atoms), (b) z(J) vs. J for two samples 
of liquid at 120K, from which glasses at (a) where created. 

A structural similarity between the glasses and liquids is 
observed at high local pressures. For lower pressures the 
environments become denser and one observes stronger 
fluctuations and lower coordination numbers. The de- 
pendence of the average radius on J is shown in Figure 
4. As a somewhat better estimate for the radius one 
can use R(J), defined as a solution of the self-consistent 
equation: 



J = z(J)U'(R(J)) + J bg 



(12) 



This radius is also given in Figure 4. Since U'(R) in- 
creases with distance, R(J) av and R(J) grow with J . 



FIG. 4. Solid line with circles R(J) av vs. J (eq. 11). 
Dashed line with squares R(J) vs. J where R(J) is calcu- 
lated from eq. 12. 

The increment in z{J) is not an obvious effect at all. 
One could argue for such effects from considerations of 
packing in the larger nearest neighbor coordination shell, 
but this will not explain the observed results. In fact 
we have obtained converse effects when other potentials 
were considered |18|, so it seems that there are various 
effects due to the potential shape. Notably the coordi- 
nation numbers are all close to 12, which is the so called 
'icosahedral' packing [4-6] . However there are very strong 
variations. One can observe the same effects in the liq- 
uid state, where the rescaled z(J) still has a square de- 
pendence on J . We note that the number of neighbors 
increases with J until it reaches a value of 13.5 and then 
one observes no more environments. The system cannot 
sustain larger z(J) and R(J) av 's because such an envi- 
ronment will collapse. 

To understand this, consider the following one dimen- 
sional example. Suppose that an atom is attached to two 
opposite walls with a potential U(R). If the walls are 
close enough there is only one stable solution with a neg- 
ative pressure. If the walls are far away there will be two 
stable solutions near the walls and an unstable solution 
at the center. The criterion for the loss of stability of 
the central solution is U"(R co ) = 0, where R co is the dis- 
tance to the wall. For the LJ potential R co — cr(26/7) 1/ ' 6 , 
whereas the distance of zero force for the single pair in- 
teraction is Rq = cr2 1//6 . Once a particle's distance ex- 
ceeds R co , it will simply attach itself to another atom. A 
glance at the distribution function P(J) shows that this 
scale is the exact cutoff for the distribution of the largest 
J's. This mechanism is probably the main mechanism 
for two-level system dynamics for such glasses. 



IV. STRESSES AND ELASTIC CONSTANTS 

Although the existence of correlations between the 
stress and the local elastic constants is expected [ |10[ , 
there is no prediction about their form. We studied the 
dependence between tensors Ai and Si and the scalars Ki 
and Ji (calculated from eqs. 0, [| || and ||, respectively) . 
We have found linear relations between Af and Sf . 
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though no such relations are found between the specific 
V and U" in eqs. and g. Those results hold despite 
strong variability between the different elements in A 
and S. A linear relation exists also between K and J: 



K t = aJi 



(14) 



Those relations are found both for glasses and for liquids. 
The small variations in the middle curve for K (see 
Figure 5) indicate that eq. [TI] is not an 'exact' mathe- 
matical identity. Numerical estimates from eq. QJ for one 
of the samples are a = —0.24 and b — 12.2. When only 
non stressed terms in eq.fflare considered, a = —0.21 and 
b = 12.0. The formal errors in the coefficients are < 1%. 
The numerical estimates from eq. QJ are a = —0.235 and 
b = 12.6. Very close numerical estimates were obtained 
for other samples that were studied (deviations < 1%). 
Results for specific eigenvalues are much more variable, 
as can be seen in Figure 5. 




(a) ^ (6) 

FIG. 5. (a) Upper, middle and lower plots are, respectively, 
the largest (IFmaxi), the average (Ki/3) and the smallest 
(Wmirii) eigenvalues of Ai vs. Ji for glass at OK. (b) Wm&Xi, 
Ki/3 and W^mim vs. Ji for liquid at 120K. 

The difference between liquids and glass is manifested 
by two main effects: the width of the plots in Figure 
5, and the variation of the definite eigenvalues of A. In 
the liquid state, fluctuations in values of the eigenvalues 
grow, so that many of the sites are unstable. 

As far as the distribution of local force constants is con- 
cerned, the K — J relation (eq. 14) provides the essence 



of the relations found in eq. (13). Therefore in the fol- 
lowing we consider only the scalars Ki and Ji. Linear 
fits for the A" vs. J are shown in Figure 6, for the same 
glasses and liquids. 
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FIG. 6. (a) K vs. J for glass at OK. Solid line - linear 
regression of the Ki — Ji curve calculated using the points in 
Figure 5. Long dashed line - linear regression of the Ki — Ji 
calculated in the same way without the stress terms in eq. 8. 
The width of the Ki-Ji plots is very narrow, with formal error 
< 1%. Dotted line - K-J plot calculated form eq. 17. (b) K 
vs. J for liquid at 120K. Solid, long dashed and dotted lines 
calculated as on (a). 

How can one explain those results? 

There are terms oc C/y in equation R that could account 
for such results, as noted in |l(J. However, such terms 
would be relevant only in systems where the pressure 
distribution is much larger (such as liquids and covalent 
glasses). In the present case U^/Rij <§C E/y, so there is 
only a small linear contribution resulting from this term. 
The K — J relation without the contribution from this 
term is shown in Figure 6. 

Since K and J are related, we can use the following 
relations to calculate them: 



J = J fp U'(r)gj(r)dr + J bg 
K = f fp U"(r)gj(r)dr + K bg 



(15) 



Expanding the potentials around some arbitrary refer- 
ence point Rq, the first order estimates for J and K are: 



J - z{J)[U'(Ro) + U"(R )SR(J)} + J bg 
K = z(J)[U"(R ) + U"'{R )5R(J)} + K bg 



(16) 



where SR(J) = R(J) av — Ro- Using eq. 16 we can esti- 
mate the K — J curve to be: 

K={J- J bg )F2(U(R Q )) + z{J)Fl(U(R )) + K bg (17) 

with 

[C/"(i?o)] 2 - U'(R )U'"(Ro) 
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A natural choice for Rq in eq. 17 is Rg = R(J) a v- The 
results for F1(U(R)) and F2(U(R)) are given in Figures 
7a and 7b, respectively. Though the corrections due to 
z(J), Fl{U(R(J) av )) and F2(U(R(J) av )) are quite large, 
the contributions of those three large dependencies cancel 
each other to create a linear curve ! The error of the 
obtained linear fit (shown in Figure 6) is less than one 
percent. 

Thus, the form of the potential is closely related to the 
structural information given by the coordination num- 
bers. Moreover, the linearity of K — J relation 
constrains the possible structure for a given po- 
tential. This holds both for glasses and for liquids (with 
our LJ potential). 





FIG. 7. (a) F1(U(R)) vs. R calculated from eq. 17 for 
a glass at OK. R is given in units of a and F1(U(R)) is in 
units of e/a. (b) F2(U(R)) vs. R calculated from the same 
equation. 

To check the extent of those results for a more gen- 
eral potential, we simulated a problem with different two- 
body potentials |18|. In preliminary simulations we have 
found that the linear laws relating K and J (as well as 
A and S) are conserved. However, the dependence of 
z( J) on J is different, it is reduced with the coordination 
number. Thus, it seems that the linear K — J depen- 
dence is a more basic law. Obviously, this hypothesis 
should be studied with a set of different potentials for 
two-body interactions, so as to get a deeper theoretical 
understanding. 



a linear relationship to the local stresses. This holds in 
glasses and is even more apparent in liquids. We have 
shown that an explanation for the linear relationship 
lies in the mutual cancellation of contributions from the 
structural arrangement of atoms and from the potential. 
Our observations indicate that the same distorted envi- 
ronments exist in the liquid and in the glass. In the liquid 
they are more distorted due to temperature, but still the 
same basic relation holds. It seems that the square de- 
pendence of the coordination number on the pressure is 
a result of this linear law. 

Our present results shed some light on the interrelation 
between local structural and elastic properties in simple 
mono-atomic glasses. In further research, a satisfactory 
theory for the dependence of the coordination number 
on J and for the linearity of the K — J curve should 
be derived using the knowledge of correlations of specific 
local pressures. 

This research is also related to some major unresolved 
issues in the physics of glasses. The obtained results for 
short range effects imply also that there will be longer 
range effects. The relevance of this implication for the 
problem of medium and long range correlations in glasses 
is discussed in a forthcoming paper [|16|. 

Knowledge of the elastic constants can provide a lot 
of information on various other aspects of glassy and liq- 
uid dynamics. In glasses, the global spectra and the de- 
pendence of the scattering on time/frequency is a major 
aspect of the analysis of the spectra. The local elastic 
constants determine the scale of the interactions in the 
global dynamical matrix, which is used to solve the global 
phonon spectra. Thus, the global dynamics is in fact de- 
termined by the local disorder. Effects observed in the 
global phonon spectra, such as the boson peak and very 
strong scattering seen in the heat conductivity, are all 
manifestation of this local disorder. 
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V. CONCLUSIONS 



In this paper we have shown how the structural and 
elastic properties of short range mono-atomic glasses are 
related. The local pressure has a nontrivial relation to 
both the radius and the coordination number. We believe 
that those relations are a result of the relation between 
the elastic constants and the local pressure. Our major 
finding is that the traces of the elastic constants have 
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